Solitonic lattices in photorefractive crystals 
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Two-dimensional spatial solitonic lattices are generated and investigated experimentally and nu- 
merically in an SBN:Ce crystal. An enhanced stability of these lattices is achieved by exploiting 
the anisotropy of coherent soliton interaction, in particular the relative phase between soliton rows. 
Manipulation of individual soliton channels is achieved by use of supplementary control beams. 
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Wide (~1 mm) Gaussian beams launched in a photore- 
fractive (PR) crystal in the self-focusing regime tend to 
break into spatially disordered arrays of filaments, owing 
to transverse modulational instabilities |l(. However, or- 
dered arrays of Gaussian beamlets (~10 /jm), launched 
in conditions appropriate to the generation of spatial 
screening solitons [2(, form much more stable solitonic 
lattices. Weakly interacting pixel-like arrangements of 
solitons that can individually be addressed are interest- 
ing for applications as self-adaptive waveguides H Q . 

Adaptive waveguides are of particular interest in all- 
optical information processing for their potential to gen- 
erate large arrays, as well as for allowing many configu- 
rations with different interconnection possibilities. Spa- 
tial optical solitons are natural candidates for such ap- 
plications, owing to their ability for self-adjustable wave- 
guiding and versatile interaction capabilities, as demon- 
strated in light-induced Y and X couplers, beam splitters, 
directional couplers and waveguides. In addition to such 
few-beam configurations, the geometries with many soli- 
tons propagating in parallel — the so-called soliton pix- 
els, arrays, or lattices — have been suggested for applica- 
tions in information processing and image reconstruction 
11110,11 Recently several groups demonstrated the 
formation of quadratic arrays of solitons in parametric 
amplifiers 8, 9] and PR media, with coherent an( i 
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FIG. 1: Experimental setup for the creation of solitonic lat- 
tices and control of individual pixels. Insets: typical input 
and output arrays. 




incoherent 4] beams. 

In this communication we combine the properties of 
spatial PR solitons to form pixel-like lattices, to in- 
vestigate experimentally and numerically the generation 
and interactions in large arrays of spatial solitons. We 
achieve improved stability of solitonic lattices by utiliz- 
ing anisotropic interaction between solitons, in particular 
the phase-dependent interaction between solitonic rows. 
We manipulate individual or pairs of solitons using inco- 
herent and phase-sensitive control beams. 

Creation of solitonic lattices requires stable non- 
interacting propagation of arrays of self- focusing beams. 
A crucial feature in the parallel propagation of PR spa- 
tial solitons is their anisotropic mutual interaction 
Because the refractive index modulation induced by a sin- 
gle soliton reaches beyond its effective waveguide, phase- 
dependent coherent as well as separation-dependent inco- 
herent interactions, such as repulsion or attraction, may 
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FIG. 2: Determining critical separation in a rectangular lat- 
tice, after 20 mm propagation, (a) Above-critical 12 x 12 
array with Aa; = 70/nn, Ay = 85/im. (b) Slightly below- 
critical 17 x 17 array with Ax — 60^im, Ay = 75/xm. (c) 
Below-critical 25 x 25 array with Ax = 50/im, Ay = 65^im. 

appear between the neighboring array elements. These 
interactions also affect the waveguiding characteristics of 
an individual solitonic channel. Therefore, the separa- 
tion between solitons and their nearest-neighbor (NN) 
arrangement, determined by the form of the lattice, must 
carefully be chosen so as to minimize all forms of inter- 
action. In our experiments, as well as numerical sim- 
ulations, we determine the critical soliton distances for 
the given crystal thickness at which the interaction be- 
comes noticeable. To achieve closer packing of solitons, 
while maintaining propagation without interaction, the 
phase relationship between different beams in the array 
is exploited. 

The creation of lattices and controlling individual pix- 
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FIG. 3: Quadratic 12 x 12 in-phase solitonic lattice after 
35 mm propagation, numerical. Electrooptic coefficient 210 
pm/V, external electric field 900 V/cm. The total size of data 
windows is 1 mm 2 and the initial FWHM size of individual 
beams 20 /im. (a) Ax = Ay = 70//m. (b) Aa: = Ay = 60/im. 
(c) Aa; = Ay = 50/im. 

els is explored in the experimental setup of Fig. ^ Laser 
beam derived from a frequency-doubled Nd:YAG laser 
emitting at 532 nm illuminates a spatial light modulator 
(SLM), which imprints the image of a spot-array onto 
the beam. The spatial light modulator in turn is imaged 
onto the front face of a photorefractive SBN60:Ce crystal 
(5 x 5 x 20 mm), which is positioned so that the propa- 
gation direction is along the 20 mm axis. To exploit the 
dominant electro-optic component w 200 pm/V of 
our SBN sample, the incident laser beam is linearly po- 
larized parallel to the c-axis of the crystal, perpendicular 
to the propagation direction. Regular patterns of up to 
25 x 25 spots, each with the diameter of about 20 /im and 
an intensity of 110 nW, are imaged onto the front face 
of the crystal. Applying an external DC electric field of 
about 1 kV/cm across the crystal and illuminating the 
crystal with uniform white light to induce an artificial 
dark conductivity, creates appropriate conditions for the 
formation of spatial screening solitons. To manipulate in- 
dividual or pairs of solitons in the lattice, an additional 
control beam is derived from the same laser. It is our ex- 
perience that large arrays of solitons can easily be formed 
in PR crystals. The larger the array the less significant 
the finite size (edge) effects. The number of pixels is 
mainly limited by the aperture of the PR crystal and 
the resolution of the inducing SLM. One can tailor the 
arrangement of beams according to specific needs. 

Propagation without mutual interaction requires suffi- 
cient initial distance between beams. It also depends on 
other parameters influencing the soliton formation and 
propagation, such as the initial beam profiles and phases, 
and the strength of the nonlinearity, but less crucially. A 
distance smaller than the critical distance for coherent 
interaction between solitons leads to their attraction and 
eventual fusion. The attractive interaction is more pro- 
nounced along the y-transverse direction, perpendicular 
to the direction of the external field. For a given propa- 
gation distance this is reflected in the deformed appear- 
ance of the lattice. Deformations may also be caused by 
the crystal inhomogeneities that can not be controlled. 
Hence the determination of the critical distance is quali- 
tative, based on the inspection of series of runs with de- 
creasing inter-soliton distances. A typical series of runs 
is depicted in Fig. [21 where the distance between soli- 
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FIG. 4: Comparing the propagation of 6 x 6 sub-critical 
square arrays with 40 fim NN separation, for different prop- 
agation distances, given in each figure. Left column an inco- 
herent lattice, middle column a coherent in-phase lattice, and 
right column a phase-engineered lattice. Other parameters as 
in Fig. H 

tons is reduced by increasing the density of packing. An 
estimate of the critical distance is found to lie between 3 
and 4 beam diameters for in-phase arrays in our experi- 
mental conditions. It is slightly larger in the y-direction 
than in the x-direction. This estimate is corroborated in 
numerical simulations. 

Our numerical simulations are based on the paraxial 
approximation to the anisotropic propagation of coherent 
optical beams in saturable PR media, 



2iknod z A 



V 2 , A 



-k 2 nQr 3 3(d x (j))A , 



(1) 



where k is the wave number in vacuum, n a w 2.35 is the 
bulk refractive index, is the transverse Laplacian, 
and A the slowly- varying envelope of the electric field of 
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FIG. 5: Propagation of square, hexagonal, and honeycomb 
solitonic lattices with 40 fim NN separation, (a)-(c) Initial 
intensity distributions, (d)-(f) In-phase lattices propagated 
for 10 mm. Note in (f) the inversion of the honeycomb lattice. 
(g)-(i) In-phase lattices propagated for 16 mm. Strong edge 
deformation is evident, (j)-(l) The same lattices, but with 
alternate rows of solitons out-of-phase, propagated for 16 mm. 
Other parameters as in Fig. [3 

the beam. The paraxial equation is augmented with the 
equation for the electrostatic potential — V<^> = E of the 
space-charge field E generated in the crystal 0, ^| 



V(£ • V ln(l + I) = E e d x ln(l + I) 



(2) 



where E e is the external field. The light intensity I is 
normalized to the background intensity. These equa- 
tions are solved together, in the manner described in Ref. 
|14|. Behavior qualitatively similar to the experimental 
was found: Parallel propagation of coherent solitons with 
varying separation can be adjusted to be almost without 
interaction, (Fig. |3| (a)) whereas for separations below 
the critical length the fusion of soliton columns is ob- 
served. (Figs. 01(b), (c)). 

Propagation of incoherent lattices was found to be 
more stable than the propagation of coherent in-phase 
lattices for identical conditions, however our intention is 
to exploit the possibility of manipulating the distribu- 
tion of phase across the lattice, to achieve closer pack- 
ing while maintaining propagation without interaction. 
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FIG. 6: Comparing numerical and experimental action of a 
control beam, placed mid-way between two pixels, (a) Initial 
configuration, (b) Numerical situation after 20 mm propaga- 
tion, (c) The corresponding experimental situation. The data 
are 70 fim pixel separation, each spot 20 fim in diameter, co- 
herent lattice, incoherent control. The control beam is three 
times more intense than any of the pixels. 
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FIG. 7: Illustrating the action of a phase-sensitive control 
beam in a coherent square lattice with 70 /im NN separation. 
Phase shift of the control beam relative to the lattice soli- 
tons is indicated on the top of each column. Numbers within 
each figure indicate the propagation distance in mm. Other 
parameters as in Fig. |3] 

The comparison between incoherent, coherent and phase- 
engineered lattice is provided in Fig. The NN sepa- 
ration is 40 /im, which is about half the critical distance 
from Fig. [21 and twice the FHWM of initial beams. It 
was found that the best results are obtained when the ad- 
jacent rows of solitons are out-of-phase. This is in agree- 
ment with theoretical expectations of the anisotropic in- 
teraction: Along rows the solitons are in-phase, and along 
columns they are out-of-phase, which makes the inter- 
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action between them universally repulsive and leads to 
enhanced stability. It is seen that the phase-engineered 
lattice in the right column of Fig. 0]is much more stable 
than the incoherent or coherent in-phase lattices in the 
left and central columns (note also different propagation 
distances). For smaller lattices, as in Fig. 0] the edge 
effects are also more pronounced. 

To further test the idea of phase engineering for en- 
hancing the stability, different lattices were formed and 
propagated (Fig. 0. It is seen that the in-phase lattices 
deform fast, and that after 10 mm of propagation, owing 
to the attractive forces along the y-direction, some rows 
are lost in the quadratic lattice, while the hexagonal and 
honeycomb lattices are strongly deformed. The honey- 
comb lattice deforms continually from the direct to the 
inverse (hexagonal) lattice and back (Figs. 03(f) and (i))- 
However, out-of-phase lattices for identical initial condi- 
tions and at a longer propagation distance of 16 mm are 
only slightly deformed. It is also evident that the edge 
effects decrease with increasing the lattice size. 

Use of solitonic arrays for applications in information 
technology requires means of manipulating individual 
waveguides, so as to combine different channels, sepa- 
rate them or induce energy exchange between the chan- 
nels. For this purpose the well-known interaction effects 
of spatial PR solitons can be exploited 0, • We uti- 
lize here a supplementary steering control beam derived 
from the Nd:YAG laser that can be focused anywhere on 
the front face of the crystal. Varying the phase of the 
control beam relative to the phase of array elements, a 
phase-sensitive coherent or incoherent interaction can be 
induced, that may lead to fusion or repulsion of different 
solitons in the array. Early experimental results on an 
incoherent control were presented in Refs. 0, 0> EH • 
Here we display the comparison with experiment, using 
an incoherent control beam acting on a pair of coherent 
pixels (Fig. |SJ). The data are as in the experiment of 
Refs. [3,ll3- Qualitative agreement is visible. 

The action of a phase-sensitive control beam is de- 
picted in Fig. [7] We chose to present the cases of and 
±7r/2 phase shifts of the control beam relative to the lat- 
tice. While for the zero phase shift the merging of the 



two channels is seen, in the case of ±7r/2 a strong energy 
exchange between the control beam and the solitons is 
observed. For tt/2 phase shift the energy flows from the 
pixels to the control beam, so that at approximately 9 
mm the control beam is maximally bright and the pixels 
are extinguished. For the phase shift of — tt/2 the flow of 
energy is from the control beam to the pixels. At approx- 
imately 10 mm the control beam is extinguished, and the 
two slightly displaced pixels are maximally bright. For 
the zero phase shift the three beams gradually merge, so 
that at approximately 15 mm only one beam is visible. 
Afterwards the cycles for the three values of the phase 
shift repeat in the inverse order, but less clearly visible, 
and some energy is transferred to the adjacent pixels (not 
shown) . 

In summary, we have investigated several aspects 
of the interaction of solitons arranged in lattices, for 
possible use in all-optical waveguides and interconnects. 
We have displayed experimentally the formation of 
large two-dimensional arrays of solitons that can be 
used to guide light at different wavelengths. Increased 
stability of solitonic lattices is achieved by the phase 
engineering of soliton rows. We have shown how to 
use a supplementary steering beam to fuse, extinguish, 
or enhance selected channels in the solitonic lattice. 
Our numerical simulations, based on the anisotropic 
model of PR solitons, are in qualitative agreement with 
experimental results. 
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